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Abstract: Objective priors for sequential experiments are considered. Com- 
mon priors, such as the Jeffreys prior and the reference prior, will typically 
depend on the stopping rule used for the sequential experiment. New expres- 
sions for reference priors are obtained in various contexts, and computational 
issues involving such priors arc considered. 
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1. Introduction 

Bayesian analysis using objective or default priors has received considerable atten- 
tion; cf. Datta and Mukerjee [17], Bernardo [6, 7] Berger and Bernardo [4], Berger 
[■)], Ghosh, Delampady and Samanta [21], and references therein. The latter book, 
in particular, contains an excellent discussion of the issues and controversies involv- 
ing objective priors, reflecting the many years of leadership of J. K. Ghosh in the 
field (along with his many coauthors). See, for example, [13, 20, 22, 23]. 
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A common objective prior is the Jeffreys prior [27], which is proportional to the 
square root of the determinant of the Fisher information matrix. The Jeffreys prior 
is quite useful for a single parameter model, but can be seriously deficient for multi- 
parameter models; this has led to preference for reference priors in multiparameter 
situations (cf. Berger and Bernardo ['>] and Bernardo [7]). 

Almost all results on objective priors have been for fixed sample size experiments. 
In practice, however, statistical experiments are often conducted sequentially, with 
a known stopping rule (cf. Siegmund [Hii] and Ghosh, Sen and Mukhopadhyay [-4]). 
Bartholomew [2] and Geisser [19] introduced the notion that objective priors for 
a sequential experiment should depend on the expected stopping time. Ye [38] 
derived the reference prior for sequential experiments when the expected stopping 
time depends on the parameter of interest only. In this paper we generalize Ye's 
result in various directions, and provide some new computational tools for use with 
priors that depend on expected stopping times. 

The paper is arranged as follows. Section 2 reviews the Fisher information 
matrix for sequential experiments with a known stopping rule, derives the Jef- 
freys/reference prior for illustrative one-parameter examples, and then provides an 
expression for multiparameter reference priors when the stopping rule satisfies a 
certain property. In Section 3, reference priors and matching priors (cf. Datta and 
Mukerjee [17]) are derived for Bar-Lev and Reiser's [I] two-parameter exponential 
family. Illustrations are given for normal distributions with several commonly used 
stopping times. 

Computation of expected stopping times is often difficult, so that utilization of 
reference priors for sequential experiments is typically challenging. In Section 4, 
an approximation to the reference prior for sequential experiments is introduced 
which is exact under some circumstances, seems to be a reasonable approximation 
in general, and allows for much simpler computation. 

2. Noninformative priors with a known stopping rule 
2.1. Notation and the Jeffreys-rule prior 

We assume that Xi,X2, . . . , is an i.i.d. sequence of random variables with den- 
sity f{x \ 6) that is regular (Walker [35]). Here is a g x 1 vector of unknown 
parameters. Let N denote a proper stopping time for the sequential experiment 
- see Govindarajulu [25] for a definition, which also is a source for the following 
well-known lemma: 

Lemma 2.1. Let 1(0) he the Fisher information matrix based on Xi. Under the 
proper stopping time N, the Fisher information based on (Xi, . . . ,Xn) is 

(2.1) r = EgiN)i{e). 

The Jeffreys-rule prior [27] for 6 is defined as the square root of the deter- 
minant of the Fisher information matrix. In the fixed sample size case, this is 
TTj{0) cx ]I(0)]^/^. For the sequential experiment, it follows from the above lemma 
that Jeffreys' prior is 

(2.2) ^}(0) (X {Ee{N)Y'^\l{e)\^'^ oc {Ee{N)Y'^T:j{e) . 

Example 2.1. Let iV^ be a random variable with a negative binomial distribution 
NB{r,p), where r is a positive integer and p e (0, 1). Let Xi, X2, ... be a sequence 
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of Bernoulli random variables with success probability p. Nr can be viewed as a 
stopping time for this Bernoulli sequence as follows: 



Nr ^ inf{n > 1 : Xi 
The probability of Nr is 



Xri 



r}. 



for fc = r, r + 1, 

An easy computation yields Ep{Nr) ~ r/p. Since the Jeffreys rule prior for a 
Bernoulli random variable is ttj{p) cx — p), it follows from (2.2) that the 

Jeffreys rule prior for the negative binomial distribution is 

r 1 
TT*j{p) cx - 7rj(p) cx 

p pV^-p 

This, of course, is well known from a direct computation with the negative binomial 
distribution, as discussed in Geisser [h\ and Bernardo and Smith ([ I !>], Example 5.14, 
p. 315). 

We next consider an example with a continuous stopping time. 

Example 2.2. Let {Z{t) : t > 0} he a Brownian motion with constant drift 9 and 
variance 1 per unit time, so Z{t) ~ N(9t, t). Let — oo < a < < 6 < oo, and let Tab 
denote the random stopping time 

(2.3) Tab = nif{t > : Z{t) < a or Z{t) > b}. 

It follows from Hall [26] that 



Ee {Tab 



(b-a) 



c'2-be 



1 



g2(6-a)e _ I 



—ab. 



if 6* 7^ 0, 
if 6* = 0. 



Note that the constant prior is the Jeffreys prior based on stopping at a fixed time 
(Poison and Roberts [29]; Sivaganesan and Lingam [31] ), from which it follows that 
the Jeffreys or reference prior for this situation is 

TT{e) = ^Ee{Tab). 

This is of additional interest because of the study in Brown [10], which showed 
that the commonly used estimate Z{T)/T, which is the posterior mean under a 
constant prior for 0, is inadmissible under estimation with squared error loss. Brown 
[10] further suggested that prior distributions which behaved like \0\~^ as \0\ oo 
were optimal for this situation. The Jeffreys/reference prior has behavior ]6']~^/^ 
as \9\ oo, and so is not of this form, but admissibility is very dependent on the 
loss function used. Indeed, it can be argued that a weighted-squared error loss is 
appropriate for this situation, and the reference prior is likely admissible for an 
appropriate weight. 
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2.2. Reference priors 

Reference priors depend on a grouping and ordering of the parameters; see Berger 
and Bernardo [4, 5]. Suppose that 6 = • ■ • j^(m)) is an m-ordered grouping, 

where the dimension of component is qi for i = 1, . . . ,m. Datta and Ghosh 
[14] considered the special case in which the (fixed sample size) Fisher information 
matrix is diagonal, with the diagonal elements being products of functions of the 
Our first result is a generalization of their result. 

Theorem 2.1. Suppose that the Fisher information matrix corresponding to a 
single observation Xi is of the form 

m m 

(2.4) 1(0) = diag(j[Gu{e^,-,), . . . , [] G™(e(.))) , 

i=l 1=1 

where Gu is a qi x qi matrix. Assume further that the expected stopping time is of 
the form 

m 

(2.5) Ee{N) = l[g,ie(,)). 

1=1 

Then the reference prior for 6 in the sequential experiment is 

ni 

(2.6) 7r^(0(i),...,0(™)) (X n[5'(^w)]''^'^«(^(i)'---'^(™))' 

1=1 

where 'n:R{6[i), . . . , 0{m)) is the reference prior based on the single observation Xi, 
given by 

m 

(2.7) 71^(0(1),..., 0(,„))=niG"(^w)r^'- 

1=1 

Proof. The proof is essentially identical to that in Datta [12], noting that, under 
(2.5), the sequential Fisher information matrix has the product structure of Datta 
and Ghosh [13, 14, 1-5]. □ 

This theorem can also be considered to be a generalization of Ye [.HS] , who con- 
sidered the case where Eg{N) depends only on the parameter of interest. 

Berger and Bernardo [•")] suggested that one should always try to use a one- 
at-a-time reference prior, where each component of the grouping of parameters 
contains only one parameter, and much of the subsequent literature has validated 
this suggestion. We thus take it as given here that a one-at-a-time reference prior is 
the desired target. The following result is an immediate corollary of Theorem 2.1. 

Corollary 2.1. Suppose that the conditions of Theorem 2.1 hold. If qi ~ 1, for 
i = l,...,77i = k, then the resulting one-at-a-time reference prior for 6 in the 
sequential experiment is 

7r^(6») oc ^Ee{N)^R{ex,...,ek). 

For later purposes, we also note another corollary of Theorem 2.1, which applies 
if the dimension of each component of the grouping of parameters has dimension 2. 

Corollary 2.2. Suppose that the conditions of Theorem 2.1 hold. If all qj = 2, 
then the reference prior for 6 in the sequential experiment is 

7r},{e) cx £;e(7V)7rfl,(0(i),...,0(,„)). 
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3. A two-parameter exponential family 
3.1. The model and reference priors 

Bar-Lev and Reiser [1] considered the following density function of the generic 
two-parameter exponential family: 

(3.1) fix 161,62) = a{x)cxp{d,Ui{x)^6iG'2{62)U2{x)-ii6i,62)}, 

where 61 < 0, fla = E{U2{X) \ (61,62)}, G^{-), (i = 1,2) arc infinitely differentiable 
functions satisfying G" > 0, and 

^{01,62) = ^6i{62G'2{62)~ 02(62)} + Gi{62). 

This is a large class of distributions, which includes, for suitable choices of Gi, G2, 
Ui and C/2, many popular statistical models such as the normal, inverse normal, 
gamma, and inverse gamma. Table 1, reproduced from Sun [ ], indicates how each 
distribution arises. 

Let Xi,X2, ... be a sequence of random variables from (3.1). The Fisher infor- 
mation per observation is 



i) 



G'{(0i) 
-6iG'^(62) 



The two parameters 61 and 62 are orthogonal in the sense of Cox and Reid [11]. 
Thus the Jeffreys prior for a single observation is 



(3.2) 7r^(0i,e2) (X VK\^G'((9i)G'^(02). 

When either 61 or 62 is the parameter of interest, it is shown in Sun and Ye [33] 
that the one-at-a-time reference priors are 



(3.3) 7:^(61,62) = ^G'{(6i)G'^(62). 

The parameter 62 is the expectation of U2(Xi). Bose and Boukai [?)] considered 
inference about 62 in sequential experimentation with the following stopping time: 

2 

(3.4) Na = inf {n > mo : r„ < nG[ (-^) }, a > 0, 

where Yn = J27=iUi(Xi) — G2{n~^J2"^^ U2(Xi)} and mo > 2 is an initial 
sample size. From Theorem 2 of Bose and Boukai [!)], we have 

TVq 1 

(3.5) lim — = — a.s. 

(3.6) lim Eg(— ^ ^ 



\0i 



Table 1 

Special cases of Bar- Lev and Reiser's [l] two parameter exponential family, where 

h{ei) = -01 + 61 iog(-0i) + iog(r(-ei)) 





Gi(0i) 


02(62) 






01 


02 




-iiog{-2ei) 




x' 


X 


-l/(2a^) 




Inverse Gaussian 


-iiog{-2ei) 


1/92 


l/x 


X 


-a/2 




Gamma 


h{ei) 


- log 9-2 


— logo; 


X 


— a 




Inverse Gamma 




- log 6*2 


log X 


l/x 


— a 
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Bar-Lev and Reiser [1 ] showed that the distribution of Yn does not depend on 
the parameter 62- So condition (2.5) satisfies when either 9i or 62 is the parameter 
of interest. The following result is immediate from Theorem 2.1 or Corollary 2.1. 

Fact 3.1. (a) The Jeffreys prior for {61,62) in model (3.1) with the stopping time 
(3.4) and when a is large is approximately 

(3.7) 7r}(0i, 02) (X ^G'l{6,)G'i{62). 

(b) The one-at-a-time reference prior for (01,6*2) in model (3.1), when either 61 
or 62 is the parameter of interest, the stopping time (3.4) is used, and a is large 
enough, is approximately 

(3.8) ^RiOiM « |^y^G'/(0i)G^'(02). 

Example 3.1. Suppose Xi, X2, . . . , are a sequence of N{^, a^) random variables. 
Then 61 = -l/2a\62 = G^^i) = -l/20i, and r„ = Y]l=l{^^ -Xn? ■ The 
stopping rule (3.4) becomes 

n 

Na = infjn > mo : n-^^{X, - X^f < nV(2a2)}. 

1=1 

So the priors (3.2), (3.3), (3.7), and (3.8) are, respectively, 

or equivalently, 

Trj{n,a) (X itr{^i,(t) oc -, 7r}(/i, cr) oc -, 7rjj(/i, cr) oc 

5.^. Probability matching priors for a sequential experiment 

Asymptotic frequcntist coverage is an often-used criterion to compare objective 
priors; see Welch and Peers [3()], Peers [28], Tibshirani [34], Datta and Ghosh [I '], 
Datta, Ghosh and Mukerjee [Ki], and Datta and Mukerjee [17] for discussion and 
references. The most common approach is to find a "matching prior," i.e., a prior 
which results in posterior one-sided credible intervals that are also accurate as 
frequentist confidence intervals. Another type of matching prior, considered by Sun 
and Ye [-V^] , is a prior such that the confidence interval based on the signed squared 
root transformation of the log-likelihood ratio is also a Bayesian credible interval. 
Almost all of the literature considers the fixed sample case for i.i.d. observations; 
exceptions are Ye [38] and Sun [32]. 

For sequential experiments involving the Bar-Lev and Reiser [1] two-parameter 
exponential family, let ln{6i,62) be the log-likelihood function of (01,6*2), given 
Xn = {Xi, . . . , Xn), and let {6„i,6„2) be the maximum likelihood estimator of 
(01,02). Write 



n n 

r„ = ^ L/i (X,) - G2 {n-i U2{X:d] . 

i=l 1=1 
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Then, on {¥„ e G[{ei)} n {n,-^ X^Li t^2(^z) e 62}, 9ni is the solution of r„ = 
Gi(^ni), and ^„2 = C^2(X0. Define 

= Gi(0i)-Gi(t^i)-G'iK)((?i-cci), wi,0ieei, 

/2('^2,e2) = G2(t^2)-G2((?2)-G^(e2)(c^2-^2), ^2,^2 6 82. 

From the convexity of Gi and G2, these two funetions are nonnegative. From Sun 
[32], the log-likelihood ratio is ^„(^„i, ^„2) - ln{0i,92) = {Z^i + ^,^2)/2, where 

Zn2) \ {^2nej2{k2,02)Y'^sgn{e2~k2) ) 

is a generalized signed square root of the log-likelihood ratio. 

Let P[0^_g.^) denote probability over Xi,X2, . . . , given ((?i,6'2), and. for a fixed 
prior 7r(6'i,6'2), let P'^{- \ X^) denote posterior probability given X„. Suppose we 
are considering a stopping time, Na, such that Na 00 almost surely as a — > 00. 
An asymptotic frequentist matching prior in this sequential setting is a prior tt such 
that 

(3.9) P^{Zn^^1 < Cl, ^Ar„,2 < C2 I XatJ 

= Piei,d2){ZNa,i < ci,ZNa,2 < C2) + 0(a~-^), 

for all ci and C2 in Pj^g^ g^)— probability. 

Suppose now that the stopping rule satisfies 

(3.10) — ^r(0), inLi. 

a 

From Sun [ ], the unique prior satisfying (3.9), and hence the unique asymptotic 
matching prior, is 



(3.11) ^U0l.02)^^r{e)G'i{e^)G'i{e2). 

As an immediate example, for the stopping time defined in (3.4), property (3.6) 
establishes that (3.10) holds; hence the reference prior given in (3.8) is also the 
asymptotic matching prior, a very desirable situation. 

Example 3.1 (continued). In deriving the sequential likelihood ratio test to see if 
{fJ.,(J^) = {fJ.0,'^0), Woodroofe [37] considered the following stopping rule, 

n 

(3.12) Na = min(b2a, inf |n > &ia : ^ - n - 7ilog(CT^) >2a|). 



where < 61 < 62 < 00 are two prespecified numbers, ~ n ^ '^^i^ii-^i ^ 
and Xn = Y^7=i -^i- Theorem 8.3 of Woodroofe [ • ] implies that 

r &2, if p'(e) < 1/62, 

p'{0), if 1/62 <P^(e)< 1/61, 
in Pj-g^ g,)— probability, as a ^ 00, where 

(3.13) p'^{9) = Gi{0i) - Gi{-0.5) - G[{-Q.5){0i + 0.5) - 0101 
= {{f^' + l)/a'+\og{a')-l}/2. 
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Thus (3.11) gives an asymptotic matching prior for this situation. Note, however, 
that the expected stopping time is not of the form (2.5), so that we cannot assert 
that this prior is also a one-at-a-time reference prior. 

4. Computation 

If Eg[N] is available in closed form, as in the examples in this paper, computation 
with any of the sequential priors can be done using common MCMC techniques. 
Hence we only consider here the case in which Eg[N] can only be computed nu- 
merically. 

4.I. Brute force computation 

All the Jeffreys, reference, and matching priors that have been discussed for a 
sequential experiment are of the form 'i>{Eg[N])Trp{6), where is some operator 
and TTi? is the corresponding prior for the fixed sample size experiment. The posterior 
distribution corresponding to this prior is 



where Xn ^ [Xi, . . . ,Xn) is the data. 

The brute force method for simulating from this posterior distribution is the 
following Metropolis algorithm: 

Step 1. Sample a proposed O' , from the fixed sample size posterior density of 0, 
which is proportional to ttp{0) YliLi fi^i I ^)- 

Step 2. Numerically estimate Eg' [N] . For instance, one could repeatedly sample N 
from its distribution given 6' , by simply repeatedly simulating the sequential 
experiment for the given 9' , observing the N that results from each simulation, 
and averaging to obtain the estimate Eg/ [N] . 

Step 3. Perform a Metropolis step: sample u uniform (0, 1) and, with denoting 
the previous value the parameter, accept 0' if 



and set 6' equal to the previous 6 otherwise. 

If one cannot directly draw from the posterior in Step 1, one could instead 
using any MCMC scheme, e.g. Gibbs sampling or Metropolis-Hastings. If doing so, 
however, be sure to iterate Step 1 many times before moving on to Step 2. This 
is because Step 2 is typically extremely expensive, as it may involve thousands of 
simulations of the entire experiment simply to compute one Metropolis acceptance 
probabihty. In situations where one dependent step is much more expensive than 
others, it pays to iterate first on the others. 



N 



(4.1) 



[OIXn) (X 'f{Eg[N])7rp{e)Y[f{X,\9) 



i=l 
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4-2. The two-dimensional case 

If using the Jeffreys prior in a two-dimensional problem or the reference prior in 
the situation of Corollary 2.2, the posterior distribution is of the form 

N 

(4.2) 7r*{9\XN) (X Eg[N]7rp{e)l[f{X,\e). 

1=1 

This allows a remarkable simplification in the computation, by introducing iV as a 
latent variable. 

To avoid confusion, we will label the latent variable as N; it is a variable with 
the same distribution as N, but is independent of N. Write the density of N given 
6 as p{N I 6). Then the joint density of {N, 9), given the data Xn = (Xi, . . . , Xjv), 
is proportional to 

N 

(4.3) p{N\e)NTTF{e)l[f{x,\e). 

1=1 

Sampling (iV, 0) from this distribution will result in 6 from (4.2), as can easily be 
seen by marginalizing over N in (4.3). 

Here is a Metropolis algorithm for sampling from (4.3). 

Step 1. Sample a proposed 0', from the fixed sample size posterior density of 0, 

which is proportional to ttp{0) YliLi fi^i I 
Step 2. Sample a proposed N' from p{N \ 9'). This can always be done by simply 

simulating the sequential experiment once, given 9' . 
Step 3. Perform a Metropolis step: sample u ~ uniform(0, 1) and, letting (TV, 0) 

denote the previous value the parameter, accept [N' , 9') if 

• ^1 

u < mm < 1, > , 

- \ iV'/ 

and set {N', 9') equal to the previous {N, 9) otherwise. (Note that, if N' < N, 
one would always accept the candidate.) 

The reason that this is a vastly more efficient algorithm than the brute force 
algorithm is that one need only simulate a single draw of N' in Step 2, whereas 
thousands of draws would be needed in Step 2 of the brute force algorithm to 
compute Eg'[N]. Again, of course. Step 1 could be replaced by any convenient 
dependent MCMC scheme. Whether one then needs to iterate Step 1 before moving 
on to Step 2 will be context dependent. 

4-3. Modified reference priors 

The most desirable prior is the one-at-a-time reference prior given in Corollary 2.1, 
resulting in the posterior distribution 

N 

(4.4) TT*{e\XN) (X ^Eg[N] nR{9)l[f{X, \ 9). 

4=1 

Unfortunately, the latent variable trick is not available for sampling from this dis- 
tribution. 
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Interestingly, however, it is frequently the case that 

(4.5) ^Eg[N]^ Eg[VN]. 

When this is the case, the latent variable trick can be applied, and the algorithm 
from Section 4.2 can be utilized by simply replacing N/N' in the Metropolis step 

with ^JWJW>. 

In the remainder of the section, wc discuss the reason that the approximation 
(4.5) often holds. The first is that the sampling distribution of N may be rather 
concentrated in a region of large N , in which case the approximation is clearly 
good. 

Example (Bar-Lev and Reiser [1]) (continued). For the stopping time Na 
defined in (3.4), it follows from (3.5) and (3.6) that 

lim Eg{y/Nja) = l/\e^\'/\ 

a — >oo 

We then have 



lim Ee\l — - lim \ Eg( — -\ 

a— too V a a^oo \ \ a / 

Example 2.1 (continued). Let Nr have the negative binomial distribution NB{r, 
p). Note that Ep{Nr) = r/p and Varp{Nr) = rp/ (1 — p)^. As r — > oo, we have 

\/ Nr/r — > in probability 

and 

Epi^KJ^) ^ ^Ep{Nr/r) = 1/vp. 

To see the difference between Ep{y^ Nr/r) and y^Ep(Nr/r) for moderate r, they 
arc plotted, as a function of p, in Figure 1 for r = 1 and r = 9. For r = 9, the 
curves are essentially indistinguishable; even for the minimal r = 1 they are quite 
close. 

It is also interesting to look at the posterior distributions for this example. In 
Figure 2, we plot the posterior densities of p for three priors 

nj{p) cx l/^/p{l-p), 
TT*jf{p) cx 1/y/p, 

(p) cx Epi,/N*J^). 

Here tta/ (p) is an approximate prior. For even the very small r = 2, the posterior 
densities under the two priors tt^ and ttm are quite close, yet substantially different 
from that under ttj. For a moderate r = 10, the posterior densities under 7r|j and 
ttm are almost identical. Note that the posterior densities of p under ttj and 7r|j are 
Beta (r, Nr — r + 0.5) and Beta (r = 0.5, Nr ~ r + 0.5), respectively. The posterior 
densities of p under ttm were computed using 5000 Metropolis samples. 

As a final indication of the similarity of the true and approximate reference 
priors in this example, and of the value of using the sequential reference priors, we 
compare the frequcntist coverage probabilities that result from their use in obtaining 



Objective sequential experimentation 



29 




Fig 1. Negative binomial example: comparison of Ep(Nr/r) and Ep(^J Nr/r) for r = 1 and 
r = 9. 



Table 2 

Coverage Probability of one-sided 5% (95%) Bayesian credible sets for the negative binomial 
Example 2.1, under the three priors ttj{p) = p(l — p), ttIjCp) = 1/(pV1 ~ 'P)> 



TTMip) = Epiy^N*/^) 



r 


P 




7r|j 




2 


0.1 


0.1142(0.9738) 


0.0516(0.9511) 


0.0487(0.9509) 


2 


0.5 


0.0002(0.9652) 


0.0010(0.9381) 


0.0008(0.9455) 


2 


0.9 


0.0001(0.9724) 


0.0003(0.9700) 


0.0000(0.9729) 


8 


0.1 


0.0751(0.9642) 


0.0474(0.9498) 


0.0465(0.9534) 


8 


0.5 


0.0552(0.9688) 


0.0522(0.9536) 


0.0568(0.9517) 


8 


0.9 


0.0000(0.9307) 


0.0001(0.9310) 


0.0002(0.9339) 


30 


0.1 


0.0617(0.9571) 


0.0508(0.9497) 


0.0516(0.9523) 


30 


0.5 


0.0556(0.9594) 


0.0512(0.9495) 


0.0525(0.9503) 


30 


0.9 


0.0426(0.9369) 


0.0438(0.9410) 


0.0442(0.9368) 



confidence intervals for p. Table 2 considers the frequentist coverage of one-sided 
5% and 95% Bayesian credible regions, based on the fixed sample size Jefi'reys' 
prior TTj, the sequential Jeffreys/reference prior 7r|j and the approximate prior ttm 
for various combination of r and p. The fixed sample size Jefi^reys' prior performs 
worse then the other two, indicating the value of using the sequential versions, while 
the reference prior and the approximate prior are almost equally good. 
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(a). Posterior Densities of p, Given (r, N_r) = (2, 5) 
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P 



(b). Posterior Densities of p, Given (r, N_r) = (10, 25) 




Fig 2. Posterior densities of p based on the priors iTj(p) = 1/ yjp{l — p), 7r|j(p) = — p), 

and 7Tm{p) = Ep{y/Nf/r) for r = 1, 10; (a) (r, A^^ = (2, 5); (b) (r, Nr) = (10, 25). 
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